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Abstract 

We study shock statistics in the scalar conservation law dtu + dxf{u) = 
0, X £ M., t > 0, with a convex flux / and spatially random initial data. We 
show that the Markov property (in a;) is preserved for a large class of ran- 
dom initial data (Markov processes with downward jumps and derivatives 
of Levy processes with downward jumps). The kinetics of shock clustering 
is then described completely by an evolution equation for the generator 
of the Markov process u{x,t), x G R. We present four distinct derivations 
for this evolution equation, and show that it takes the form of a Lax pair. 
The Lax equation admits a spectral parameter as in [35], and has remark- 
able exact solutions for Burgers equation (/(u) = ii^/2). This suggests 
the kinetic equations of shock clustering are completely integrable. 
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1 Introduction 

We consider the scalar conservation law 

dtu + d^f{u) ^0, xeM,t>0, u{x,0) =uo{x), (1) 

with a strictly convex, flux / and initial data uq that is a stochastic process 
in X. The basic model of this type was introduced by Burgers in his study of 
turbulence. He considered f{u) = and white noise initial data [14, 44, 46]. 
While this model fails to describe turbulence in incompressible fluids, it still 
serves as a widely useful benchmark for theoretical methods and computations 
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in turbulence. It also has fascinating links with combinatorics, mathematical 
physics and statistics, some of which we describe below. 

Our main contribution in this article is to develop a consistent kinetic the- 
ory that describes completely the clustering of shocks and the statistics of the 
random process u{x,t) for convex /, and initial data uq that are Markov pro- 
cesses in X with only downward jumps. Our approach includes exact solutions to 
Burgers turbulence as an important special case. In addition to being compre- 
hensive, our approach clarifies the essential features of the problem, and reveals 
the role of a Lax pair that describes the shock statistics. In order to describe 
these results and place them in the context of past work, let us first describe 
the structure of solutions to (1) for deterministic uq, and then some important 
exact solutions to Burgers turbulence that motivated our work. 



1.1 The Hopf-Lax formula 

Let us first recall the notion of the entropy solution to (1) [20, 33]. Charac- 
teristics for (1) are lines in space-time along which u is constant. If a unique 
characteristic connects (a,0) to {x,t) then u{x,t) — uo{a). However, there may 
be many characteristics that pass through {x,t). Well-posedness is resolved by 
adding a small dissipative term euxx to the right-hand side of (1) and pass- 
ing to the limit £4-0. This was first carried out by Hopf in his pioneering 
work on Burgers equation [27], and later generalized to convex / by Lax [32]. 
Let /*(s) = sup„g]u (us — /(w)) denote the Legendre transform of / and call 
Uo{s) = uo{r)dr the initial potential. We define the Hopf-Lax functional 

I{s;x,t) = Uo{s) + tf* (^^y (2) 

The 'correct' characteristic through (a;, t) is given by the variational principle 

a(x, t) = sup < s e K : I{s; x, t) = inf /(r; x,t)>, (3) 

reK J 

We will always assume that Uq has no upward jumps and / satisfies the growth 
condition 

lim I{s; x,t) = oo. (4) 

This ensures the infimum of 7 is a minimum, and that a{x, t) is finite. We write 

a{x, t) = arg+minseR/(.s; .t, t). (5) 

This is the Hopf-Lax formula for the inverse Lagrangian function a{x,t). The + 
in (5) denotes that we choose a{x, t) to be the largest location if I is minimized 
at more than one point. Of particular importance is Burgers equation with 
f{u) = u'^/2. In this case, (5) is called the Cole-Hopf formula, and takes the 
form 

/ \ + • ( TT I \ (.'£■ — .s)^ 1 , , x — (i(.r.t) 
a{x,t) = aig'^mmseR<Uo{s) -\ — — j- , u[x,t) = (6) 
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For fixed t, a defined by (5) is non-decreasing in x. If a; is a point of continuity 
of a(-,t), the velocity field is given implicitly by 

f'{u{x,t)) = ''-f^'\ (7) 

u{x, t) is well-defined because /' is continuous and strictly increasing. In par- 
ticular, if a{x, t) is constant for x in an interval, we obtain a rarefaction wave, 
a may also have upward jumps. These arise if the minimiun of I is attained 
at more than one point. Jumps in a give rise to shocks in u. The left and 
right limits u± = u{x±,t) exist and the velocity of the shock is given by the 
Rankine-Hugoniot condition 

u(x,t) = ^ _ =:[/]«_, u+- (8) 

We stress that the entropy solution ^(a;, i) has only downward jumps in x. This 
will play a key role in our analysis. 



1.2 Exact solutions for Burgers turbulence 

Let us now suppose uo is a stochastic process in x. The solution develops 
shocks that are separated by rarefaction waves (see for example, the compu- 
tations in [40]). The shocks move with speeds given by the Rankine-Hugoniot 
condition, and cluster when they meet. The Lax equation we derive describes 
this process. In order to explain this, we first focus on two exact solutions for 
Burgers turbulence. 



Levy process initial data on a half-line 

We assume f{u) = and 

/ X / 0, a; < 0, ... 

where X is a Levy processes with only downward jumps (a spectrally negative 
Levy process) . A particularly interesting case is when X is a standard Brownian 
motion (that is, a Brownian motion with E(X^) = x). This problem was solved 
formally by Carraro and Duchon [15, 16] and rigorously by Bertoin [9]. The key 
to their solution is a closure property of (6): if uo{x), x > is a Levy processes 
with only downward jumps, then so is u{x,t) — w,(0,f), x > 0. These processes 
are characterized by their Laplace exponent ip{q,t) defined by 

^ (^ei(u{x,t)-u{o,t))'-j ^ ^x^{q,t) ^ x,9,t>0. (10) 

Thus, the problem is reduced to determining the evolution of tp. Remarkably, 
ip{q, t) satisfies Burgers equation in the new variables q, t\ 

5tV + Va,V = 0, V(9,0) = Vo(9). (11) 
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If X is a Brownian motion, ■^0(9) = 9^, and we obtain the self-similar solution 

i){q,t) = ^xP,{qt), V'*(9)= 9+^-^9+ J- (12) 

Various explicit formulas are summarized in [38]. 

The interpretation of (11) in terms of shock clustering is as follows, tf) satisfies 
the celebrated Levy-Khintchine formula 

/•oo 

i}{q,t)= / (e-9«-l + gs) At(rfs), g>0,i>0, (13) 

where At(ds) is the jump measure at time t. The evolution of ip by (11) also 
induces evolution of the jump measure A^. But the jmnps in the process u(-, t) 
are the shocks, which evolve in a simple manner: shocks move at speed given 
by the Rankine-Hugoniot condition (8) and stick, conserving momentum, when 
they meet. Equation (11) captures this process. It is equivalent to the fact 
that At satisfies a kinetic equation, Smoluchowski's coagulation equation with 
additive kernel, which describes the evolution and coalescence of shocks. A 
derivation of the kinetic equation from this perspective may be found in [38]. 
An excellent survey of several links between stochastic coalescence and Burgers 
turbulence is [11]. 



White noise initial data 

Here we must characterize the law of u{x, t) when the initial velocity is white 
noise. Precisely, let us suppose that the initial potential Uo{x) = aB^ where 
B is a standard two-sided Brownian motion pinned at the origin and cr a fixed 
scale parameter. This problem also arises in statistics, and it was in this context 
that Groeneboom first characterized the law of the process u{x,t), x [26]. 
He showed that for every t > 0, u{x,t), x e R is a Markov process with only 
downward jumps (a spectrally negative Markov process) that is stationary in x. 
He then computed the generator of this Markov process explicitly in terms of 
Airy functions. Here is a brief summary of his solution. 

By (6) and the scaling invariance of Brownian motion, we see that 

aa{x,t) = {at)^ai yx{at)~'s ,1) , Ua{x,t) = a'^t~'Sui lx{at)~^ ,1) , (14) 

where = denotes equality in law and the subscript a refers to the variance of Uq. 
It is simplest to state the formulas under the assumption that = 1/2. Let 
denote the process M2-i/2(x, 1). The generator of is an operator defined by 
its action on a test function ip in its domain as follows: 

^^(y)=lim^?^M^^ll^^, (15) 

xiO X 

where Ej^ denotes the law of the process with uq = y. Groeneboom showed that 
A is an integro-differential operator of the form 

ry 

Aip{y)=ip'{v)+ \ {ip{z)-if{y))n^{y,z)dz. (16) 
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The jump density n* of the integral operator is given explicitly as follows. 



n,{y,z)=^^K{y-z), y > z, (17) 

where J and K arc positive functions defined on the line and positive half-line 
respectively, whose Laplace transforms 

/OO /^OO 
e-''yj{y)dy, k{q) = e-'^yK{y)dy, (18) 
-QC Jo 

are meromorphic functions on C given by 
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J(9) = TT73' HQ) = -'^:j-2^ogAi{q). (19) 



Ai denotes the first Airy fimction as defined in [2, 10.4]. Our normalization 
of J and K differs from [26] in two aspects. First, the above definition of K 
is suggestive of Dyson's formula in the theory of inverse scattering [34, p. 273]. 
Second, the choice = 1/2 helps avoid several factors of 2^/^ while stating the 
main formulas. 

We shall return to this sohition at several points in this article. The one- 
point and two-point distribution functions can be computed once the generator 
is known. For example, p{y)dy = P{ux € {y,y + dy)) = J{y)J{—y)dy. The 
distribution of the shock sizes and the two-point distribution is given in [23]. 
Tauberian arguments yield precise asymptotics of these distributions. 



1.3 Lcix equations 

Despite the elegance of the solution procedure for Levy process data, it does not 
apply to Burgers equation with broader classes of random initial data (e.g. white 
noise), or to the general scalar conservation law (1). Our goal in this article is to 
develop a kinetic theory that describes shock clustering for (1). Our work builds 
on the links between stochastic coalescence and Burgers turbulence [11, 38], 
exact solutions to Burgers equations with white noise initial data [23, 26], and 
recent work on kinetic equations for Burgers turbulence [17]. We amplify these 
remarks briefly. 

Groeneboom used the variational principle (6) at t = 1/2 to compute the 
generator as in (16)-(18). Girsanov's theorem is crucial in the analysis. Most 
past work in the turbulence community is also based on a similar point of 
view. This begins with Burgers' analysis which was eventually completed by 
Frachebourg and Martin [8, 14, 23, 30]. However, the focus on white noise 
initial data and the use of Girsanov's theorem obscures an understanding of the 
dynamic process of shock clustering for general initial data (see for example, 
the concluding remarks in [10]). 

The idea of seeking kinetic equations for shock evolution can be found in pio- 
neering early work of Burgers [12, 13], but this approach remained undeveloped 
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for several deeades. More recently, E and Vanden-Eijnden introduced a hierar- 
chy of master equations for the statistics of velocities and velocity gradients in 
Burgers turbulence [22]. However, these equations are not closed, and a sophis- 
ticated dimension reduction is needed to extract scaling exponents from these 
equations. Frachebourg et al showed that the hierarchy of n-point functions 
for ballistic aggregation can be closed at the level of the 2-point function [24]. 
Most recently, Chabanol and Duchon showed formally that if a statistical so- 
lution to Burgers equation preserves the Markov property, then one can derive 
evolution equations for the generator of the Markov process [17]. Moreover, 
they showed that Groeneboom's solution yields a self-similar solution to this 
evolution equation. 

This is the starting point for our work. Our viewpoint is as follows: rather 
than seek a specific exact solution to Burgers equation with specific initial data, 
we look for a class of natural stochastic processes whose structure in x is pre- 
served by the entropy solution to the general scalar conservation law (1). Here 
it is the class of spectrally negative Markov processes (in x). This reduces the 
study of shock statistics to an evolution equation for the generator of the Markov 
process. Our main results are: 

• a closure theorem: Suppose / is strictly convex and C^. If uq is a spec- 
trally negative Markov process in x, so is the entropy solution to (1) (see 
Theorems 2 and 3) . This provides some rigorous justification for our ap- 
proach. 

• new kinetic equations: We derive kinetic equations that describe the evo- 
lution of the generator of the Markov process for arbitrary convex /. The 
kinetic equations are equivalent to a Lax pair. Moreover, the Lax equa- 
tions admit a spectral parameter as in Manakov's integration of the Euler 
equations for a spinning top in [35]. This provides strong evidence 
that the kinetic equations are completely integrable. 

• consistent derivations: We derive the Lax equation from four different 
perspectives. Aside from being a stringent test on the consistency of the 
kinetic equations, this also provides a unified treatment of disparate meth- 
ods in the literature on statistical hydrodynamics. 

In order to explain the Lax equation, let us assume that u{x,t) is a station- 
ary, spectrally negative Feller process in x whose sample paths have bounded 
variation. As in the Levy-Khintchine formula, such a process is characterized 
by its generator A{t), which acts on test functions ip e C^(]R) via 

Afiy) = b{y, t)^'{y) + f {^{z) - ip{y)) n{y, dz, t). (20) 

For fixed t > 0, b{-,t) £ C(M) and n{y, •, t) is a measure on (— oo, y) that satisfies 
/(I A |j/ — z\'^)n{y,dz,t) < oo for each y e R. These terms correspond to the 
drift and jumps of the process respectively. There is no diffusion term because 
sample paths are solutions to a conservation law and necessarily have bounded 



6 



variation [20] . This expression for ^ is a special case of the general form of the 
generator of a Markov process, as given in Theorem 3.5.3 in [7] (see also [41], 
Eq. (2.1.13)). 

We now introduce an operator B associated to A and the flux function /. 
This operator is defined by its action on test functions as follows: 

B^iy) = -ny)b{y, t)^'{y) - f ^^^^ ^^'^ {^{z) - ^{y)) n{y, dz, t). (21) 

One of our main results is that the evolution of A is given by the Lax equation 

dtA=[A,B] = AB-BA. (22) 

An evolution equation for the generator was first derived in [17] for the flux 
f{u) = However, the equations in [17] were not written in the simple 

form above, and it was not clear if such equations could even be derived for 
arbitrary fluxes /. We show that while f{u) = u^/2 is certainly special, much 
of the structure of the problem relies only on the convexity of /. 

We present four different derivations of (22): (1) as a compatibility condition 
for martingales in x and t; (2) from kinetic theory as in [38]; (3) using BV 
calculus as in [22]; (4) using Hopf's functional calculus as in [17]. The shortest 
(and most heuristic) of these is the first, and goes as follows. 

The main observation is that we have a two-parameter random process, and 
formally B may be viewed as the 'generator' in t. To see this, fix x <E R and 
consider the random process u{x,t), t > 0. Then the multipliers ~f'{y) and 
— [f]y,z in (21) simply correspond to the evolution (in t) of the drift and jumps. 
Indeed, if the path u{x, t) is differentiable at x with u{x, t) = y, dxU = b, then 
dtu = —f'{y)b by (1). Similarly, if we have a shock connecting left and right 
states y and z, then the shock speed is given by (8). 

This may be understood more precisely using Ito's formula for jump pro- 
cesses [7]. For fixed t, the random process u{x,t) satisfies the stochastic diff'er- 
ential equation 

du{x,t) = b{u{x-))dx + / {z — u{x-))N{u{x-),dz,dx), (23) 
Js. 

where the Poisson random measure N{y,dz,dx) has intensity n{y,dz)dx. We 
write the conservation law as du(x,t)dx + df{u{x,t))dt = 0, and apply Ito's 
formula to obtain 

du{x,t) = -f{u{x-))b{u{x-))dt- / [f]u(^x_)^^{z- u{x-))N{u{x-),dz,dt). 

Jm 

This shows immediately that B should have drift —f'{y)b{y) and jump measure 

^[I\y.zn{y. dz). This holds rigorously when / is decreasing. In this case, all 
shocks move to the left, u{x-,t) = u{x,t-), and u is also Markov in t, with 
generator B given by (21). 
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If A and B arc generators, then the one-point distributions satisfy Kol- 
mogorov's forward equations 

d^p = A^p, and dtP^B^'p. (24) 

We now seek martingales in x and t. To this end, fix (xo,to)j and consider the 
processes (p {u{xo + s, to)) and (p {u{xo, to + s)) with s > 0. These processes are 
formaUy martingales if ip solves Kolmogorov's backward equations 

d^if + Aip^O, dt(p + Bip = 0, (25) 

in the domain {x,t) € [a;o,oo) x [to,oo), and y,z S R. If the compatibility 
condition pxt = Vtx holds for a sufficiently rich class of functions p, we obtain 
the general Lax equation 

dtA-d,B = [A,B]. (26) 

If the process is stationary in x, d^B vanishes and we obtain (22). In the form 
(26), the Lax equation is akin to zero curvature conditions in integrable systems. 



1.4 Kinetic theory 

When we expand the commutator in (22), and separate the evolution of the drift 
b and the jump measure n we obtain a kinetic equation that describes shock 
clustering. The drift satisfies the differential equation 

dtbiy,t) = -f"iy)b\y,t). (27) 

Note that the drift does not depend on the jump measure. The jump density 
n{y, z, t) dz = n{y, dz, t) satisfies the kinetic equation 

dtn{y, z, t) + dy {nVy{y, z, t)) + d, {nV, {y, z, t)) (28) 
= Q{n, n) + n (([/],,, - /'(y)) dyb - bf"{y)) . 

Here the velocities Vy and Vz in (28) are given by 

Vy{y, z, t) = {[f]y,z - f'{y)) b{y, t), V^iy, z, t) = ([/],,, - f{z)) b{z, t), (29) 

and the collision kernel Q is 

Q{n,n){y,z,t) = [ {[f]y^^ - [f]w^z)n{y,w,t)n{w, z,t) dw 



- [ {[f]v,z - [f]z,w) n{y, z, t)n{z, w, t) dw 

J —oo 

/V 
{[f\v,w ~ [f\v,z) n{y, z, t)n{y, w, t) dw. (30) 
-oo 

We have assumed for convenience that the jump measure has a density, but the 
above equations extend naturally to general jump measures. 
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In the next section we derive these equations from the perspective of kinetic 
theory. We consider single shocks and rarefaction waves as building blocks, 
and use this to derive a natural Boltzmann-like equation. This will yield the 
equations above. We then show that these kinetic equations are equivalent to 
the Lax equation (22). That calculation reflects the fact that operators of the 
form (20) formally constitute a Lie algebra. 

1.5 The broader context of our work 

We conclude this introduction by connecting our work with some other prob- 
lems in mathematical physics and statistics. We pay particular attention to 
connections with integrable systems. 

1.5.1 The spectral curve and complete integrability 

We first point out the role of a spectral parameter in analogy with Manakov's 
treatment of the Euler equations for geodesic flow on so{n) with a left-invariant 
metric [35] . Let M and J\f denote multiplication operators acting on the domain 
of A, defined by 

M^{y) = yip{y), Afip{y) = f{yMy). (31) 

It is clear that M and JV are diagonal operators. We now use the deflnitions 
(20), (21) and (31) to find 

[A,Af\-[M,B]=0. (32) 

This observation allows us to introduce a spectral parameter e C in the Lax 
equation. We use (22) and (32) to obtain 

dt{A- nM) = [A- tJ.M,B + ijJ\f], iJ.eC. (33) 

If A, B were n x n matrices, it would follow that the spectral curve (Riemann 
surface) 

r = {(A, ^) e \dct{A - Aid - iiM) =0}, (34) 

is fixed by the evolution. In Manakov's work, this is the crucial observation that 
yields the existence of additional integrals for Euler 's equations in so{n), n > 4. 
These integrals are simply the coefficients of the characteristic polynomial above. 

More broadly, the observation that (22) admits a spectral parameter reveals a 
close relation with a large class of completely integrable systems (including KdV, 
the Toda lattice, geodesic flows on so{n) and ellipsoids, and the integrable PDEs 
of random matrix theory) . The complete integrability of all these flows may be 
obtained in a unifled way via a general splitting theorem for Lie algebras [3]. 
This connection also sets the stage for the application of powerful methods from 
algebraic geometry to integrate (26) explicitly for every convex / [4]. We will 
address this in later work. 
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1.5.2 Burgers turbulence and random matrices 

The solution to Burgers turbulence with spectrally negative Levy process data 
(see § 1.2) is obtained from (22) as follows. Suppose u{x,t) — u{0,t), a; > is 
a spectrally negative Levy process with bounded variation and mean zero. If 
we denote the jump measure At, then b{t) = Jg°° sAt{ds), and the generators A 
and B take the form 

POO 

A{tMy) = / {ip{y-s)-ip{y) + sip'{y))At{ds), (35) 
Jo 

B{tMy) = -yA{tMy) + ^j^ {ip{y - s) - ^{y)) sAtids) . 
In particular, when ip{y) = e^^, Ke{q) > 0, we have 

^(f)e«« =^(g,^)e«^ B{t)e^y ^ - ( y^{q,t) + ^-dgA e'^y . (36) 



We substitute (36) in (22) to obtain (11). 

We next note that the solution (12) can be mapped to Wigner's semicircle 
law in the theory of random matrices. Dyson observed that the eigenvalues of a 
standard matrix valued Brownian motion Mt in the group of n x n symmetric, 
hermitian or symplectic matrices satisfy the stochastic differential equation [21] 



where 1 < k < n are independent Brownian motions, and ,6 = 1, 2 or 4 for 
the ensembles above. That is, the eigenvalues behave like repulsive unit charges 
on the line perturbed by independent white noise. The law of large numbers 
for this ensemble is as follows. As n — >■ oo the spectral measure of n~^/'^Mt 
converges to Wigner's semicircle law: 



fjit{dx) = -^^^4* - a;2 dx, \x\<2Vt. (38) 

Moreover, the Cauchy transform of /xt, 

g{z,t) = j j^pLt{dx), 2; e C\[-2\/i,2\/i], (39) 

solves Burgers equation with a simple pole as initial data [5, §4.3.2]. That is, 

dtg + gd,g = Q, g{z,Q) = -. (40) 

z 

More precisely, g solves (40) in the slit plane €\[—2^/t,2\/t], and has the form 
5(z,t) = i=5*(^), g.{z) = \(z~^^4), \z\>2. (41) 
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It now transpires that the self-similar solution (12) can be transformed to (41) 
by a simple change of variables. 



= 9*{z), z = 2 + - 



or 



(42) 



In his beautiful thesis [29] , Kerov found a deeper interpretation of (40) based 
on the representation theory of the symmetric group. He introduced a Markov 
process for the growth of Young diagrams {Plancherel growth) ^ and derived (40) 
in a mean-field limit. More general initial conditions may also be included in 
(40) and he showed that the evolution of g by Burgers equation is equivalent 
to a kinetic equation for /xt [29, Ch. 4.5]. Thus, the transformation (42) links 
/it to At, and Plancherel growth with Smoluchowski's coagulation equation (see 
(13)). We do not have a deeper (i.e. stochastic process) explanation of this 
relation yet. 

Airy functions and the Painleve transcendents arise in the scaling limit of 
fluctuations from Wigner's law at the edge of the spectrum. The fluctuations 
are given by the celebrated Tracy- Widom distributions involving a solution to 
Painleve-II [42]. We now point out that the function I = j' / j {j as in (18)) 
solves the Riccati equation dl/dq = —q + P and is therefore an Airy solution 
to Painleve-II [1, Ch. 7]. This is used to verify that Groeneboom's solution 
satisfies the kinetic equation (22) in §6. 

1.5.3 Shell models of turbulence and their continuum limits 

It is also of interest to consider (1) on the half-line a; > with random forcing 
at a; = 0. This problem arises as a continuum limit of shell models of turbu- 
lence [36, 37]. Shell models are lattice equations of the form c'fe = Jk-i — Jk, 
k = 1,2,.... Here Cfc > models the energy in the A;-th Fourier mode (shell), 
and Jfc the flux from shell k to fc + 1. is expressed in terms of and its 
nearest neighbors by a suitable constitutive relation. The main question is to 
understand how randomness spreads through the system under the assumption 
that Jo is a prescribed random forcing (a stationary Feller process, for example). 
In particular, it is of interest to understand whether the system has a unique 
invariant measure. In the continuum limit, these questions may be treated by 
our approach. We note that the derivation of (26) does not depend on assump- 
tions of stationarity and is independent of boundary conditions. Therefore, (26) 
holds in the domain x,t > and must be augmented with a boundary condition 
on the line x = 0. In particular, as t ^ oo, invariant measures are solutions to 



with a boundary condition at a; = matching B and the generator of the forcing. 

In contrast to a time-correlated boundary forcing, it is also possible to con- 
sider a white-in-time forcing in the bulk (a Feller process in x and ^-correlated 
in t) which is independent of the initial data. For Burgers' equation, this is 



(43) 
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a particular case of forced Burgers turbulence. The generator T of the forcing 
simply appears as an additional term in the Lax equation 

dtA-d^B=[A,B\ + T. (44) 

An interesting exact solution to (44) for Burgers equation forced with a two- 
sided Brownian motion in space has recently been obtained in [18] . 

1.5.4 Applications to statistics 

In statistics, (6) with Uq a two-sided Brownian motion first arose in the follow- 
ing estimation problem [19]. Suppose Xi, . . . , X„ are independent, identically 
distributed (iid) samples from a distribution with a smooth unimodal density p 
with mode m and finite variance. Let us consider a naive 'binning' strategy to 
estimate the mode m. We fix a bin width w and count the number of samples 
Nn{s) = #{Xk S {s — w,s + w)} in the bin centered at s. We estimate the mode 
by TO„ = arg+maxsAf„(s). Chernoff observed that the fluctuations m„ — m are 
0(n^^/^). Precisely, for suitable c{p,w) > 0, the rescaled random variables 
c'n}/^{mn — m) converge in law to Chernoff 's distribution 

Z = arg+maxs {Uo{s) - s^] . (45) 

The quadratic term — arises from the Taylor expansion of p around m - the 
first order term vanishes since m is the maximum of p. By the symmetry of 
Brownian motion, it is clear that Z has the same law as a(0, 1/2). 

This is not an isolated result: 'cube-root' fluctuations appear naturally in a 
wide class of estimation problems [6]. Kim and Pollard proved functional limit 
theorems for several such estimators, with the law of the limit characterized by 
arg+max^gj^d {Uq{s) — jsj^} where Uq{8) is a continuous Gaussian process in 
pinned at the origin [31]. These correspond to solutions to Burgers equation in 

with random initial data Uo, but this connection has not been explored in 
the literature. 

1.6 Outlook 

Let us conclude by pointing out some significant shortcomings in our work. 
Much of this article relies on formal calculations. But these calculations are 
often interesting, and it seems more fruitful to present them in a transparent 
and suggestive manner, rather than as rigorous statements burdened by techni- 
calities. The most significant gap is that we do not prove (22). This is because 
our closure theorem is not strong enough. We only establish that the entropy 
condition preserves the Markov property. In order to rigorously establish (22) 
it is necessary to prove that the entropy condition preserves Feller processes. 
This is an assertion of regularity, whereas we only establish measurability. This 
is closely tied to establishing a satisfactory well-posedness theory for (22), and 
will be addressed in forthcoming work. A suitable well-posedness theorem would 
also yield a probabilistic proof of the existence of a two-parameter family of self- 
similar solutions to (22). These solutions are generated by considering the flxix 
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functions f{u) — \u\p /p, I < p < oo and an initial potential that is an a-stable 
spectrally negative Levy process. It seems challenging to prove this analytically 
starting with (22). 

Also, while OUT work yields a deeper imderstanding of Groeneboom's solu- 
tion, it does not, as yet, constitute an independent proof of his results. We have 
only been able to verify that his solution satisfies (22). We have been unable 
to derive it using (22) alone. We hope to address this in future work using 
techniques from integrable systems. 

The rest of this article is organized as follows. We first derive (27) and (28) 
from the standpoint of kinetic theory in the next section. This is followed by the 
rigorous closure theorems. We then derive the Lax equations by BV calculus 
and by Hopf's method. Finally, we consider Groeneboom's self-similar solution 
for Burgers equation in Section 6. 

2 Kinetic equations 

2.1 Introduction 

In this section, we assume the velocity it is a stationary Feller process in x, and 
derive the kinetic equations of Section 1.4. We use the evolution of a single shock 
and rarefaction wave to derive a Boltzmann-like equation for the evolution of the 
density of shocks. We conclude this section by showing that the Lax equation 
(22) is equivalent to the kinetic equations (28)-(30). The main observation is 
that the space of operators of the form (20) is formally a Lie algebra. 

2.2 Conservation of total number density 

Let p{y,t) denote the stationary 1-point density, i.e. p{y,t)dy = P{u{x,t) G 
{y,y + dy)), and F{y,z,t) denote the total number density, i.e. the expected 
number of jumps per unit length from states y to z. Then 

F{y,z,t)=p{y,t}n{y,z,t). (46) 

The total number density changes because of a flux of shocks and shock colli- 
sions, and we have the general conservation law of Boltzmann-type 

dtF + dy (FVy) + (FV,) = C{F, F). (47) 

Here Vy and denote 'velocities' in the {y, z) 'phase space', and C is a binary 
collision kernel. This is the general structure of the equation. We now derive 

the evolution equation for 6, the velocities Vy, Vz, and the collision kernel C, 
based on elementary solutions to the scalar conservation law (1). 

2.3 Decay of the drift 

First consider how affine data evolves under (1). Let u{x,t) solve (1) with 
uo{x) = ao+ l3ox. For x,t 0, we have to leading order u{x, t) w a{t) + I3{t)x, 
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so that 

dtu^a + ^x, d,f{u)^ f{a)P + f"{a)P^x. (48) 
We now balance terms in (1) to obtain 

a = ~f'{a)f3, p^-f"{a)p\ (49) 

The second equation expresses the decay of rarefaction waves when /3o > 0. The 
connection between this elementary solution and the generator of the process is 

the following. The drift coefficient h{y, t) corresponds locally to an affine profile 
as above with a = y and (3 = b. Thus, the second equation above is simply (27). 

2.4 Decay of shocks and Vy,Vz 

The 'velocities' Vy and Vz in (y, z) space arise because of the decay of shocks. 
In order to derive these velocities, we fix U- > u+ and consider piecewise affine 
initial data 

Let s{t) denote the path of the shock. Then by (8), for small t 

S{t)^[f]u_,u+{t + 0{t)). (51) 

s is also given by the kinematic condition 

s{t)=a±{t) + f{u±)t + o{t). (52) 

where a± denotes the left and right inverse Lagrangian points at time t. Thus, 

a±{t) = ([/] f'{u±))t + o{t). (53) 

Since ?i_ > 7i+ and / is convex, we see that — f'{u^) < 0. Similarly, 

[/]«_,«+ — f'iu+) > 0. As a consequence, a shock initially connecting states u± 
decays to a shock connecting states 

^± + ([/] f'{u±))b±t + o{t). (54) 

This decay gives rise to a flux of F. To flrst approximation, the flux is linear in 
F, and of the form (47) with the drift velocities given by (29). 

2.5 The collision kernel C{F, F) 

Binary collisions of shocks occur at a rate determined by the Markov property (in 
x) and the relative velocity of shocks given by (8). In order to simplify notation, 
we suppress the t-dcpendence for the process u and write for w,(.x, t). We also 
denote a shock connecting states u- and u+ by {u-,u+}. Shock clustering 



involves the following events: 

growth : {y, w} + {w, z] = {y, z], z <w <y, (55) 

decay: {y, z} + {z,w} — {y,w}, —oo<w<z, (56) 

decay: {Wjy} + {y, z} = {w, z}, y < w < oo. (57) 
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The computation of rates for these events is similar. To be concrete, let us first 
consider (55). Fix z < w < y and consider small Aa;i > 0, Ax2 > 0. Then 
formally, by the Markov property 

dP{uo = y, UAxi = z) f=i p{y)n{y, z) dy dz Axi, (58) 

and similarly, 

dP {Uo = y, UAxi = W, UAa;i+Aa;2 = z) (59) 

« p{y)n{y, w)n{w, z) dy dz dw Axi IS.X2- 

The relative velocity of these shocks is [f]y.w — [f]w.z to leading order. We thus 
set = {[f]y,w — [f]w,z) At to compute the number of collisions in time At. 
We now sum over w in the range z < w < y to obtain the growth term in 
CiF,Fy. 

Ci := / P{y)n{y,w)n{w,z){[f]y^y,-[f]w^;,) dw. (60) 

J z 

The computation for the events (56) is similar. We now find the decay terms 
C2 := - / p{y)n{y,z)n{z,w) ([/]j^,^ - [f]z,w) dw, (61) 

J —00 

and 

/■oo 

C3 := - / p{w)n{w,y)n{y,z) {[f]w,y - [f]y,z) dw. (62) 
Jy 

2.6 Kinetic equations for n 

We have now defined all the terms in (47). In order to obtain an equation in 
terms of n alone we use Kolmogorov's forward equations to eliminate the 1-point 
distribution p{y, t). The first equation in (24) is now = A'^p since the process 
is stationary in x. The second equation in (24) implies 

P{y)dtn{y, z) = dtF{y, z) - n{y, z)B^p{y), (63) 

where 

S^p(t/) = dy {f {y)h{y)p{y)) - / {p{w)n{w, y) - p{y)n{y, w)) [f]y,yj dw. (64) 

Jr 

(It is convenient to denote the domain of integration by M, noting that it is 
actually a half-line because n{y, 2) = for y < z.) 

To isolate the main cancellations on the right hand side of (63), we note that 
the integral term in C3 — n{y, z)B^p{y) is 

Hv^z) ([f]y,z / p{w)n{w,y)dw -p{y) / n{y,w)[f]y,^dw] . (65) 
\ Jr Jr / 



15 



The first integral above can be simplified further. Since A'^p = 0, we also have 

/ p{w)n{w,y)dw = p{y) / n{y , w) dw + dy {b{y)p{y)) . (66) 
Therefore, we may rewrite the expression in (65) as 

p{y)n{y,z) ( {[f]v,z-[f]v,w)n{y,w)dw + n{y,z)[f]y,^dy{h{y)p(:y)) . (67) 

J —oo 

We now collect all terms on the right hand side of (63) using (47), (29), (60), 
(61), and (67). Wc then have 

P{y)dtn{y, z) = p{y)Q{n, n){y, z) (68) 
+niy, z) ([/],,.9, {h{y)p(y)) - dy (f {y)b{y)p{y))) (69) 
-dy {F{y, z)Vy{y, z)) - {F{y, z)V,{y, zj) , (70) 

where Q{n,n) denotes the collision kernel in (30). Finally, we use (46) and (29) 
to obtain the sum of (69) and (70): 

= Piy) in {[.f]y.z - f'{y)) dyb - bf"{y)) - dy {nVy{y, z)) - 8, (nV^Xy, z)) . (71) 

Under the assumption that the Markov process u has a strictly positive station- 
ary density, we may cancel p{y) on both sides of (67). We then have the kinetic 
equations (28). 



2.7 Equivalence of the LcLx equations and kinetic equa- 
tions 

The equivalence of the Lax equations (22) and the kinetic equations (28)-(30) 
follows from the algebraic structure of operators of the form (20) and (21). The 
space of operators of the form 

A^{y) = b{y)ip{y) + / n{y, z) {<p{z) - ip{y)) dz, (72) 
Jr 

with 6, n smooth, n satisfying appropriate integrability conditions, and the 
bracket [•, ■] is formally a Lie algebra. That is, if i = 1, 2, 3 are operators as 
above, then [Ai , A2] is an operator of the same form, and the following Jacobi 
identity holds: 

[[Ai,A2],A3] + [[A2,A3],Ai] + [[A3,Ai],A2]=0. (73) 

Wc do not assume that b and n are positive or that their siipport is a half-line. 
Both assertions rely on tedious, but direct calculations. We will omit the proof 
of (73), and simply summarize the calculation for [^1,^2]- We have 

[Ai,A2My) = {bib'^ - b[b2) <p{y) + / + <j) (^(2) - <^{y)) dz, (74) 

Jr 
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with 



bi{y)dyn2 - b2{y)dyni + {bi{z)n2 - 62(^)^1) , 



(75) 



and 




+ / ni{y,z)n2{y,w) - n2{y,z)ni{y,w)dw. 



+ / n2{y, z)ni{z,w) — ni{y, z)n2{z,w) dw 
Jr 



I {m{y,w)n2{w,z) - n2{y,w)m{w,z)) dw (76) 

R 



We apply (74) to (22) as follows. Let A\= A and A2 = B, with A and B as in 
(20) and (21). We then find immediately from (75) that the drift coefficient of 
[A,B] is 

-b{y) {f{y)b{y))' + f {y)b{y)b\y) = -ny)b\y) (77) 

as in (27). Similarly, we use (75) to obtain the second, third and last term in 
(28). Finally, we substitute the jump measures in (76) to obtain the collision 
kernel Q{n,n). 

3 Closure theorems 

3.1 Introduction 

In this section wc show that the entropy solution to (1) preserves the class of 
spectrally negative Markov processes. The Markov property of u was implicitly 
used by Burgers, and made explicit in [8]. Our work is based on previous results 
by Bertoin [9] and Winkel [45] . We extend these results to general convex fluxes 
/ and a large class of noise initial data. 

3.2 Splitting times 

The main technical tool we need is the decomposition of Markov processes at 
random times so that the past and future are conditionally independent. This 
certainly holds at a stopping time. However, it also holds for a broader class of 
splitting times. We use the following theorem on the preservation of the Markov 

property at a last-passage time. 

Theorem 1 (Getoor [25]). Consider a cddldg strong Markov process Xg. Let 
M CR be a fixed set and let L = sup{s e M : e M} be the end of M. Then 
the post-L process {Xs}s>l is independent of {Xs}s<l given Xl- 

Furthermore, the transition semigroups of the pre- and post-L processes can 
be explicitly determined from that of X. It is an important and subtle fact 
that these transition semigroups are different from that of the original process. 
In particular, the semigroup of the post-L process is the same as that of X 
conditioned to never hit M again. 
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3.3 Closure for spectrally negative initial velocity 

Wc work within the canonical framework. Let uq = {n,T,Ts,UQ{s), (Iq) be a 
spectraUy negative Markov process on the space f2 of cadlag paths on (—00,00) 
endowed with the Skorohod topology. Here, /xq is the law of the coordinate 
process uq{uj: s) = cj(s) on and = (j(uo(^)i r < s) is the natural filtration of 
Mo extended to be right-continuous and complete. Note that in our work x plays 
the role of "time" as traditionally used in the theory of stochastic processes, t 
acts as a parameter in the following discussion, and will be dropped from the 
notation when convenient. 

We rewrite the Hopf-Lax functional and inverse Lagrangian function as 

I{s;x,t) = 1^ (^uoir) - {fT' (^)) (78) 

a{x,t) = a.Tg^ mml{s;x,t). (79) 

To ensure that (79) is well-defined and finite, we assume (4) holds a.s. This 
is certainly true if Wq is a stationary random process and /* grows fast enough 
at infinity. In what follows we also make use of the fact that for fixed t, a{x,t) is 
an increasing function in x. Since u{x,t) and a{x,t) are related by (7), in order 
to obtain the closure property we only need to show that a{x,t) is a Markov 
process. We prove the following: 

Theorem 2. Let he a spectrally negative strong Markov process such that 
(4) holds a.s. Under the law /iq and for any fixed t > 0, the inverse Lagrangian 
process a{x, t) is a Markov process. Thus, u{x, t) is a spectrally negative Markov 
process. 

Proof. Fix X e M. Without loss of generality it suffices to take t = 1. We also 
suppress the dependence on t in the notation for clarity. The proof consists of 
three steps, the first two of which are entirely deterministic. 

1. Dependence structure of {a{x))xeM.' Since a{x) is increasing, for /i > 

a{x + h) ^ arg+ min /(s; x + h) 

= arg+ min {/(a(x); x + h) + (/(s; x + h) — I{a{x); x + h))} 

s>a{x) 

/ {uo{r)-{f')-'^{x-r))dr} . (80) 

Ja{x) J 

(a(x + h) — a{x))h>o therefore depends on uq only through {u(){s))g^^(^^y The 
same argument shows that {a{x — h) — a{x))h>Q depends only on (Mo(s))s<a(a;)- 

2. Downward jumps determine UQ{a{x)): If uq has only downward jumps 
then dsl{a{x);x) = 0. That is, 

uo{a) = {fr\x-a{x)). (81) 
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This may be seen easily by sketching a picture. Upward jumps in uq give rise to 
a potential with a corner that is convex, and downward jumps give a potential 
with a corner that is concave. In the second case, the minimum of I can never 
be achieved at the corner. Thus, it is always obtained at a point of continuity, 
implying (81). 

3. a{x) is a splitting time for uq: It is not obvious that a{x) is a Markov 
time for uq. It is not a stopping time since the event {a{x) < s} is not Ts- 
measurable. Indeed, a{x) is the last time (78) is minimized and depends on 
{uo{r))r>s as well. a{x) is, however, a splitting time for a certain functional of 
Uq. To see this, first define 

m(s; x) = min/(r; a;), D{s; x) = {I — m){s; x). 

r<s 

For simplicity, denote d{s;x) = {f')~^{x — s). Dropping the exphcit dependence 
on X from the notation, consider the process {uo,d,D). Since D{s) is entirely 
dependent on (uo(r))r<s, {uo,d,D){s) is J^g-measurable. For any J^-stopping 
time ^ 

m(C + s)=m(OAf min /(r)) = 7(0 + |-7?(0 A min (7(^ + r) - 7(0)| 

yj<r<{+s J 0<r<s J 

and 

m +s) = im m) - i-m) a min + r) - m)] . 

0<r<s J 

For s > the increment 7(^ + s) — !{£,) depends on only through wo(0 by 
the strong Markov property of uq. Therefore, (mq, d, D) is also a cadlag, strong 

Markov process. 

We now show that a{x) is a splitting time for {uQ,d,D). By definition 
(79) and step 2, a{x) is the last time s that {uo,d,D){s) hits the fixed set 
{{y, y,0) : y gM.}. Thus, we can use Theorem 1 to split uq at a{x) into 

{uo,d,D)s^a{x) and {uo,d, D)s>aix)- 

As a consequence, 

{uo{s))s<aix) and iuo{s)) s>a{x) 

are conditionally independent given d{a{x)) — that is, since d is invertible, given 
a{x). 

This shows that a is a Markov process (note that the law of the increments of 
the process may vary with x). Therefore, u{x) = {f')~^{x — a{x)) is a spectrally 
negative Markov process. □ 

Remark 1. Theorem 2 reduces to Bertoin's closure theorem [9, Thm. 2] if 
f{u) = u^/2 (Burgers flux) and uq is a spectrally negative Levy process. To see 
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this, use (81) in (80): 



/•a{x)+s 

/ {uo{r) - Uo{a{x)) + r - a{x)) i 

J a(x) 



a{x + h) — a{x) = arg+ min I I ("o('') — Uo{a{x)) + r — a{x)) dr 

s>0 I ./^(^) 



= arg"'"min<^ / {uo{a{x) ^ r-) - uo{a{x)) + r) dr \ . 

Since uq has independent inerements, the integrand above is independent of 
(uo(s))s<a(x) has the same law as UQ{r) — mo(0) + r. So, the law of a{x + 
h) — a{x) is independent of a{x — h) — a{x))h>o and does not vary with x. Note 
that it is necessary to use f{u) = v? /2 to obtain this result, and to show that 
the increments of a{x) are identical in law to those of the first hitting process 
T{x) = inf{s > : uo{s) + s>x}. 



3.4 Closure for noise initial data 

We now prove the analogue of Theorem 2 for noise initial data. Recall that the 
Hopf-Lax functional and inverse Lagrangian function satisfy 

I{s;x,t) = Uo{s) + r(^^ (82) 

a(x, t) = arg^ min/(s; x, f). (83) 
We prescribe the law of the potential as follows. Let 

III («^) 

where y, Y are independent copies of a spectrally negative additive process 
starting at and X- denotes the limit from the left. We denote by Mq the law 
of Uq on the space of cadlag paths. If X has stationary increments it is a Levy 
process, but it is not necessary to assume this. In addition, assume the growth 
condition (4) holds almost surely. 

Theorem 3. Suppose Uq = {^l,G,Qs,Uo{s),Mo) is a two-sided spectrally nega- 
tive process with independent increments and satisfies (4) a.s. Then under Mq 
and for all t > 0, u{x, t) is a spectrally negative Markov process. 

Proof. Again, we fix x e M, take t = 1, and drop the f-dependence in our 

notation. 

1. Dependence structure of {a{x))xeR'- As before, for h > 
aix + h) = arg"*" min J(s; x + h) 

= a{x) + arg"'" min {I{a{x) + s;x + h) — I{a{x); x + h)} 
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Let Xg denote the strong Markov proeess I{s:x + h) and let a = a{x). The 
increment a{x + h) — a{x) therefore depends on Uq only through the increment 

2. Spectral negativity of X: Since Uq is spectrally negative, so is the process 
Xg. Explicitly, 

Xs_ A Xg = Xg. 

3. Markov property of X at a: The inverse Lagrangian a is the last time 
Xs- A Xg hits its ultimate minimum m = Xa- Using spectral negativity, we 
need only consider Xg. Let Mg = mmr<g Xr- Then a{x) is the last time s that 
{XgjMg) hits the fixed set {{y,y) : y € K}. Theorem 1 then implies {Xs)s>a is 
independent of {Xs)s^a given m. 

We now make use of the particular form of the semigroup of the post-a 
process, for which we refer to [25, Thm. 2.12] (see also Theorem 5.1 in [39] and 
the subsequent remark). Precisely, given a and m, the law of {Xs)sya is identical 
to that of an independent copy of Xg started at a with Xa_ = m and conditioned 
to live above m. By the invariance of Uq under translations in state space (due 
to the property of independent increments) this implies that {Xs)s>a is identical 
in law to a copy of Xg started at a with value Xa = and conditioned to live 
above 0. To summarize, although the post-a process no longer has independent 
increments, the increment Xa+g — Xa is still independent of m given a. 

By steps 1 and 3, the increment a{x + h) — a{x) is independent of {a{x — /i) — 
a{x))h>o given a and m. Since a{x + h) — a{x) only depends on X through the 
increment Xa+s — Xa, given a it is independent of m as well. Therefore, a{x) 
is a Markov process and u{x) is a spectrally negative Markov process. Notice 
that if Uo is not spectrally negative or does not have independent increments, 
then a{x + h) — a{x) depends on Uo{a-) or Uo{a) in addition to a, destroying 
the Markov property. □ 



4 BV calculus and the Lax equations 

In this section, we derive the Lax equation (22) using the conservation law (1) 
and the Vol'pert chain rule for BV functions. This calculation is similar in 
spirit to [22]. We assume that for every t > 0, u{x,t) is a spectrally negative, 
stationary Feller process. The generator of u is given by (20) and the operator B 
is defined by (21). In addition, we define one and two-point operators as follows. 
We associate a linear functional V to the stationary one-point distribution p{y, t) 
via 

Tiip) = E {ip{u{x, t))) = f ^{y)p{y) dy. (85) 

(Here and in what follows, it is convenient to suppress t in the notation.) Sim- 
ilarly, given h > we denote the transition kernel for the process u by qn and 
define the associated transition operator 

{Qh^M = [ qhiy,zMz)dz. (86) 
Jr 
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4.1 The 1-point function 

Wc recall that the entropy solutions to (1) are in BVioc (R+ x R) . Vol'pert showed 
that one may extend the chain rule to such functions in a natural manner. Let 
u± = u{x±,t) denote the right and left limits of u{x, t) and for any test function 
ip, consider the composition (fi{u{x,t)) and set 

[ip]= I ip{u.+P{u+-u.))dp. (87) 
Jo 

Then for every smooth test function with compact support, the entropy solution 
to (1) satisfies [43, p.248] 

dt(fi{u{x,t)) = -[ip][f]dxu{x,t), x€m,t>0. (88) 

The law of the 1-point function is determined by E {ip{u{x, t)) for arbitrary (p. 
We use equation (88) to obtain 

dtE{ip{u{x,t))) = -E{[ip][f]dMx,t)). (89) 

The left hand side of (97) is {dtP)ip. We must determine the right hand side. 

To this end, fix .t G R and for /i > let us denote u^h = u{x — h,t). Our 
calculation relies on the following unjustified interchange of limits: 

Emmu{x,t)) = hmE fl i^iu+) - viu.h)) -/(^-^)) ') . (90) 

hiO \h U+ — U-h J 

We assume (90), and show that 

hm E (\ (^(n+) - ^{u.h)) (/(^+)- /("-^)) ') = -V (B^) , (91) 
As a consequence, 

{dtV)ip = VB^, (92) 

which is the formal forward equation (24). 

Wc now establish (91) for a monomial f{u) = u'". Since B depends linearly 
on /, (91) then also holds for polynomials and by approximation for fluxes 
with polynomial growth. For f{u) = u™ we expand the jump term 

f{u+) - f{u.,) ^ u\'-^ui. (93) 

Therefore, the limit in (91) is a sum of terms of the form 
{ip{u+) - ip{u-h)) u"!^^~^u\ = 
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The limit of each of these terms is 



lim - ,p{u-h)) u"l^^- 



= V (y^-i-^" (^(t/V) - [y^])) . (94) 

The term A{y^ip) — ^pA (y^) may be computed using (20). First the drift term 
is 

ym-i-j ^^^yj^y _ ^^^yjy^ ^ 6^™- V'. (95) 
Similarly, the jump term simplifies to 

ym-l-j f ^^y^ ^^^J (^(^^ _ ^(y)) (96) 



We now sum over all terms in the expansion (93) to obtain 
hm E - (v'(w+) - ^(u-h)) 



h^O \h U+ — U-h 

= J2'P{y"'-'-'{Ay''p)-'pAy')) 

m — 1 

~^z^ {ip{z) - ip{y)) dz 



(m-l 



= V f6my"-V' + / n{y,z)^^^^-^ {ip{z) - ^{y)) dz) = -VBip. 

\ Jr z-y ) 

4.2 The 2-point function 

Fix X e M and a > 0. Let ^ and i/' be two test functions. The law of the 2-point 
function is described completely by 

E {lp (u(x, t)) i) (u{x + a, t))) = V (ifQcip) . 

The BV chain rule extends to a product rule as follows 



where ^ 

^ix,t) = -{ip{x,t_) + ip{x,t+)). 

We combine the BV chain rule and the conservation law to obtain 

dt(p{u{x,t)) = -[(p][f]dxu{x,t), dti){u{x + a,t)) = -[tlj][f]da:u{x + a,t). 
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Therefore, 

dt {^{u{x, t)i,{x + a, t)) = -^[^'] [fPMx, t) - [f]dMx + a,t). (97) 

We compute the expected value of each of these terms in turn. First, as earher, 
the main assumption is 

E(^[^'][/']a,«(.x,t)) (98) 

h^o h \ u+ — u-h \ ^ 

where tp-^- = tp{u{x + a, t)) and ip-h = ip{'U'{x + a — h, t)). 

We can do away with the mean value. To be explicit, we write the above 
expectation as 

P{y)qh{y,z) I — ^— ^ — {ip{z) - ip{y)) I dzdy. (99) 

We see that the last term converges to Qa'4'{y) as /i — ^ 0. Therefore, 

E{i]y'][f']dMx,t)) 

= lim U /' /("+) -/(»-^) _ Q^^{u+) 

h^O n \ — U-h 

= lim 1e /" /("+) -/(»-^) (<^Q„^(w+) - ^Q„^{u_h)) (100) 
h^o n \ u+ — u- 

- y(tx-,,) -^^"+^~^^"-^^ {Qai^{u+) - Qai^{u-h))) , (101) 

= -r (B i^Qa^) - ^BQ^ij) , (102) 

where we used (91) to compute the limits in (100) and (101). 

We now compute the second term on the right hand side of (97). As in the 
calculation above, we can do away with the mean value, and we have 

E(^[/][/']9,u(x + a,0) 

= lim^ / / / p{y)^{y)q„_^{y,z)^^^^ ( fM ~ fi^) (^(^) _ ^(^))') dwdzdy. 



''^OJeJrJr h \ w — z 

As in the computation of (91) we find that the limit of the innermost integral 
is Bip{z). Therefore, 

Eilp[iP'][f']dMx + a,t)) = -V{ipQaB^). (103) 

We combine (97), (102) and (103) to find 

dt {V (^QaV')) = V{B (^Q«^) - ^fiQ«V + ■ (104) 

But the left hand side is simply 

{dtV)'pQci^ + V{^dtQoci^), (105) 
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and by (92) 

(dtVWQo^-ij^ = VB ((^QaV) • (106) 
We now combine (104), (105) and (106) to obtain 

V{ipdtQo,i>) = V{^{Qo,B-BQo,)i>). (107) 

This equation holds for all <p and ij) in the domain of A. Thus, we may write 

dtQa = {QaB-BQ^) = [Q^,B], (108) 

and in the limit a \.Q we have 

dtA = [A,B]. (109) 

5 Hopf 's method 

In this section we derive the Lax equation (22) following Hopf's method [28]. 
This method was used by Chabanol and Duchon to derive their kinetic equa- 
tions. We show that the calculation works for any convex flux /. 

Though our calculations are mainly formal, we begin with the canonical 
framework of § 3. Recall that the initial data wq has law on O. Let u{t; uq) : 
[0,oo) — >■ be a weak solution of (1) with initial data uq, which induces a 
sequence of probability measures {nt)t>o on fJ. To derive an equation for the 
flow of measures [iJLt)t>o, make the assumption that u is differentiable in t so 
that the conservation law can be written for all ip G (R) as 

dt{u{t; uo), ^) - {f{u{t; uo)), <f') = 0. (HO) 

Here, ' = dx and (•,•) is the standard duality pairing. Define the Hopf charac- 
teristic functional jit of the law fit 

At(^) = / e^<"'^Vt W (111) 

which evolves according to 

JQ. JQ 
JQ 

i{fiu),^')e'^^''^>lit{du). (112) 

Any set of probability measures {fJ,t)t>o on U that satisfies (112) for all G 
(M) is defined to be a statistical solution of the scalar conservation law — in 
particular, the flow of measures generated by the entropy solution is a statistical 
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solution. Assuming that all moments of fit are finite, the exponential can be 
expanded in series (denoting (ix„ = UJ^^ciXj) as 



oo ^„ n 

12 ~ '[lu{Xj)ip{Xj)dXn- 

n=0 -^K" j=l 



Let 



E 



'/it 



i=l 



Y\_'Pi{u{xi))iJ.t{du). 



Substituting the expansion for the exponential into (112) yields the infinite 
hierarchy 



n=0 -^K" 



'[[u{Xj)(f{Xj) 



n=0 



E 



f{u)(p'{x)J[u{Xj)(p{Xj) 



dxdxn- (113) 



We simplify this hierarchy as follows. If g : M" — >■ M satisfies g{xi,. . . , Xn) = 
9{Xcr{i), ■ ■ ■ Xcr(n)) for all permutations cr of {1, . . . , n}, then 

/ g{xi,...,Xn)d-Kn = n\ g{xi,...,Xn)dKn- (114) 

•^R" Jxi<X2<---<Xn 



n+1 



We choose 

g{xi,...,Xn+i) = -^E^, 

k=l 

in (114) and substitute in (113) to obtain 



f{u{xk))(p'{xk) Yl 'u{xj)v>{xj) 



CO „ 

n=0 Jxi<---<x„ 



Y[u{Xj)ip{Xj) 



3 = 1 

n+1 



dXr, 



E„ 



= E^""7 E 

Re-indexing the l.h.s. of (115) and subtracting the r.h.s. gives 



/(w(Xfe))(^'(Xfe) ]J u{Xj)(fi{Xj) 



dx„. (115) 



n=l •^xi<-<a;„ [^^-^^ 



fe=i 



/(w(a;fe)) Jl u{xj) 



iitk J 

(116) 
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Finally, the statistical hierarchy (116) can be considerably simplified with 
the closure assumption that {nt)t>o is the law of a stationary Feller process with 
the one and two-point operators V and {Qh)h>o defined in (85) and (86). With 
hi = Xi- Xi-i and Qi = Qhi for i = 2, . . . , n, 



E 



JJ(y9i(w(a;i)) 



p{dui,t)ipi{ui) / qh2iui,dU2,t)ip2iu2) ■ ■ ■ / qh^{Un-l, dUn,t)iPn{Un) 
= 'P(piQ2(p2 ■ ■ ■ Qnfn- 

Assume the transition measures qh are differentiable in h. Transferring the 
derivative on the test function ip'{xk) onto qhk{uk-i,duk) in (115) and using 
integration by parts, the statistical hierarchy is 



ra=l 



J|m(.Xj) 



fe=i 



Yl^PiXjJdXn = 0. 



Note that the boundary terms in the previous equation vanish due to cancella- 
tion of terms and the fact that G C^{M.) has compact support. The density 
of tensor products of the form (p{xi, . . . , x„) = 11^=1 fi^j) in the space of test 
functions on {x = {x\, . . . , Xn) G M" : a;i < • • • < Xn} implies that for n S N 
and Xi < X2 < ■ ■ ■ < Xn, the following infinite set of equations hold: 



n "(^j' 



fe=i 



f{'>^{Xk))J{u{Xj) 



(117) 



In terms of the operators = dh^Qi (denoting hi = Xi, Qi = V, Qn+i = I) 
this is 



^{PuQau • • • Qkf{u)Q'k+iU ■■■QnU 
fe=i 

- VuQ2U ■ ■ ■ Q'kf{u)Qk+lU ■ ■ ■ QnU}, 

or equivalently, in terms of the generator A (using Q- = AQi) 



VuQ2U ■■■Qk [f{u),A] Qk+iu ■■■Qr, 



k=l 
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One obtains the evolution equation for the 1-point function by taking hi — > 
for all i (i.e., Qi = I for all i), or the equation for the 2-point function by taking 
hi+i = a for Z e {1, . . . , n — 1} and /ij ^ for alH ^ Z + 1 (i.e., Qi+i = Qa and 
Qi = I for alliy^l + 1): 

71 

dtE^, [(«(xi))"] = [f{u),A]u^->' (118) 

fe=i 



k=l 

n 

+ E (119) 

k=l+l 

Using the expression (20) for the generator A, (118) and (119) simplify: 

^*^y" = -E^y'"'(W(2/)y""'+ / n{y,dz){f{z)- f{y))z^-^\ 
k=i ^ J* J 

= -r\b{y)f'{y){y-y+ [ n{y,dz)l^^tJM(^,n _ 

I Jm. z-y ) 

dtiVy'Qy''-') 

= -E^y'"' +_^ri(2/,rf^)(/(^) - f{y))z'-'Qa.z^-'^ 

E / n{y,dz){f{z)-f{y))z^->'\ 

r 111 V •'R J 



k=l+l 



-'p{b{y)ny){y'yQ^y^-'+ [ n{y,dz)l^^^^^{z' - y')Q„z^ 
I Jm. z-y 



-Vy'Q^{b{y)f{y){y--'y + / n(y, dz)^^^^^^(^"-' - y""'). 
I z-y 

In terms of the operator B defined in (21) the above expressions read 

dtVy^ = VBy", dt{Vy'Q^y^-') = V{B{y'Q^y^-')-y'BQ^y^-'+y'Q^By^-'). 

Since B is a linear operator, we approximate ip,tp & C^{R) by polynomials to 
obtain the evolution equation for the 1- and 2-point operators: 

dtVif = VBip, dtiVipQaij) = V{B{ipQ„ilj) - ipBQai^ + fQaBi^). 

Therefore we have the evolution dtA = [A, B] exactly as in § 4. 



28 



6 Groeneboom's solution 

In this section we verify that the generator 

A(tMy) = \^\y) + J' -l^n^yf/^ zt'/^) {^{z) - ^{y)) dz (120) 

satisfies the Lax equation (22) when /(«) = and n* is given by (17)-(19). 

In this case, the evolution of the drift (27) is simply b = —6^. It is clear that 
b{y,t) = is a solution. The equation for the jump density now takes the 
form ^ 

dtn{y, ^' *) - ^(2/ - ^) i^y'^ - ^zfi) = Q{n, n) (121) 



with the collision kernel 



Q{n,n){y, z,t) = — - — / n{y,w,t)n{w, z,t) dw (122) 



/ \ y ~ w , , , 
-n{y,z,t) / n{z, w, t) dw 

J -oo ^ 

/y iij — ^ 
——n{y, w, t)n{y, w, t) dw. 
-oo ^ 



We substitute the ansatz 

u = yt^/^^ v = zt^/^, s = u-v, n{y,z,t) =t-^/^n^{u,v), (123) 
in (121) and collect terms. We must then verify that 

3 + r 3)j{v)^\3 V J{u) 3 if(s) 

-is) - s——{v) + '-r^{v)) + V^(w) - s——{u) 



K ' ' V J J 7 V J 

Here x plays the role of a dummy variable in the convolutions J * (xK). Ex- 
plicitly, J * {xK)(u) = Jjjj J{u — x)xK{x) dx. Note also that the terms in (124) 
are evaluated at the arguments u, v or s = u — v respectively. The functions J 
and K are related by the following identities: 

x'^J = K*J + J' (125) 

^J*{xK) = x{K*J) + J (126) 
x^K = ?,x{K*K)+AxK' + 2K. (127) 

(The argument of each function in these identities is x). We substitute these 
identities in (124), and collect the three terms corresponding to the arguments 
u^v and s and find that (remarkably!) each of them reduces to a polynomial, 
and the sum of these three polynomials vanishes. 
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The identities (125) (127) are proven using the Laplace transform, (18), and 
the definition Ai"{q) = qAi{q). First, we observe that / = j'/j = — Ai'/Ai solves 
the Riccati equation 

l' = -q + p, (128) 
As a consequence of the definition of in (18) we have 

/' = k/2. (129) 

Therefore, we also have 

k' = 21" = -2 + All' = -2{1- Ik). (130) 

Thus, I and k solve an autonomous system. Equation (128) may be rewritten 
in terms of j and k as 

j" -{q + k)j = 0. (131) 

This is equivalent to (125). Next, if we rewrite (130) in terms of j and k, we 
find immediately that 

\3k' = {jk)' -j. (132) 

This is equivalent to (126). Finally, we differentiate (130) twice and use (128) 
to eliminate I to find 

k" = Alii -Ik) (133) 

k'" = 3(A;2)' + 4gfc' + 2A:, (134) 

which is equivalent to (127). 

Here is the Painleve property: differentiate (128) to find = 21^ — 2lq — 1. 
Now rescale r = —q2^/^, w{t) = 2^^/'^l{q) to see that w solves the second 
Painleve equation with parameter 1/2 [1]. 

l^ = 2w^ + WT+l. (135) 
ar^ 2 
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